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ABSTRACT 

In the context of the VLA-COSMOS Deep project additional VLA A array observations at 1.4 GHz 
were obtained for the central degree of the COSMOS field and combined with the existing data from 
the VLA-COSMOS Large project. A newly constructed Deep mosaic with a resolution of 2.5" was 
used to search for sources down to 4cr with la w 12 ^Jy/beam in the central 50'x50'. This new catalog 
is combined with the catalog from the Large project (obtained at 1.5" x 1.4" resolution) to construct a 
new Joint catalog. All sources listed in the new Joint catalog have peak flux densities of >ha at 1.5" 
and/or 2.5" resolution to account for the fact that a significant fraction of sources at these low flux 
levels are expected to be slighty resolved at 1.5" resolution. All properties listed in the Joint catalog 
such as peak flux density, integrated flux density and source size are determined in the 2.5" resolution 
Deep image. In addition, the Joint catalog contains 43 newly identified multi-component sources. 
Subject headings: cosmology: observations - radio continuum: galaxies - surveys 



1. INTRODUCTION 

In recent years, several cosmological deep fields have 
been im aged at 20cm (e. g., R ichards 2000; Bondi ct ^ 
200a ICondon et all [20031 IHopkins et al.l 12005 



Sevmour et all 120041: iNqrriset a.1.1 120051: iHuvnh et all 



2005t iFomalont et al.l 120061: iSimpson et al.l 120061: 



lyison et all 120071: iSchinnc rer et al.l 120071: iMiller et all 
2008; Ow en fc Morr ison 2003) providing a few thousand 
radio sources down to flux limits of a few 10/iJy. 
These deep radio imaging data are sensitive enough to 
detect star forming galaxies with star formation rates 
of several 10 to 100 M0yr~^ out to and beyond a 
redshift of z~l. Similarly, radio galaxies can be seen 
out to redshifts of z '^ 5 and the most luminous ones 
even well into the epoch of reionization. Thus, deep 
radio images in conjunction with deep imaging data 
at X-ray, optical and infrared wavelengths are ideal 
to investigate the dust-unbiased star formation, the 
evolution of radio(-loud) AGN, as well as the population 
mix of radio sources in the first place. 

In order to study the cosmological evolution of galax- 
ies and black holes, it is not only important to over- 
come the effect of cosmic variance (e.g. by studying 
a large enough area) but also to understand the ef- 
fect of large scale structure on the evolution (e.g. by 
covering a large contiguous area). To address the sec- 
ond effect in particular, the Cosmic Evolution Survey 
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(COSMOS)^ collaboration has conducted panchromatic 
imaging and spectroscopy of an equato rial field with 
a size of 2 deg^ (for an overview, see iScoville et al.l 
I2007bf) rangin g from X-ray XMM-Newto n and Chan- 
dra (iHa singer c t^l2007l: lElvis et al.ll200l . UV GALEX 
(jZamo iski ct al. 2007*), optic al and near-in frared ground- 
based (iTaniguchi et al. 20071: iCapak et alli2007'). o ptical 



HST (IScoville et al.ll2007al:lKoekemoer et al.ll2007D . mid- 

to far-infrared Spitzer (iSanders et al.ll2007|) . millimeter 

Bertoldi et al.ll2007t IScott et al.l 120081) and radio VLA 



Schinnerer et al.l 120041 . 120071 ) imaging to extensive op- 



tical spectroscopy using the VLT/VIMOS and Magel- 
lan/I MACS instruments (jLillv et al.l 120071 : iTrump et al.l 
l2007f ). Most of these datasets are now publicly available 
from the COSMOS archive at IPAC/IRSA^. 

The VLA-COSMOS survey at 20 cm is part of the 
overall imaging effort and its scientific goals and moti- 
vation have been described in detail by iSchinnerer et aLJ 
(|2007f ). Initial observations from a pilot project test- 
ing the mosaicking strat egy and giving a first source 
catalog are presented by ISchinnerer et al.l (|2004f l. As 
a large NRAO/VLA program, the VLA was used in A 
and C configuration to cover the entire COSMOS field 
resulting in an image with uniform noise properties in 
the central Ixl d eg^ a nd an average rms of 10.5 /iJy. 
ISchinnerer et al.l (|2007f ) provide a detailed description 
of the survey set-up, the data reduction, as well as 
the testing and construction of the final VLA-COSMOS 
Large p roject catalog (here after: Large catalog). Subse- 
quently, iBondi et al.l (|2008D derived the completeness of 
the Large catalog and also analyzed the effect of band- 
width smearing on the derived source flux densities to 
obtain the source counts. Although the VLA-COSMOS 
Large Project dataset has been used for several scien- 
tific results on, e.g. the faint radio population, the 
radio-derived star formation rate density, the radio AGN 
opulation an d stac king o f high-z galaxy popu lations 
Smolcicet allllM l2009allbl : ICarilh et al.ll2007l l2008h . 
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the need for deeper radio imaging data became appar- 
ent during the search for radio counterparts to milhme- 
ter sources from the COSMOS MAMBO mapping data 
(jBertoldi et al.|[2007| ): the Large project only provided 
counterparts for about half of the mm-sourccs. Thus, 
the Deep project was initiated with the aim of doubling 
the integration time for the central seven pointings which 
fully cover the MAMBO 1.2 mm map of the COSMOS 
field. These new observations are described here. 

The paper is organized as follows: After a description 
of the new observations and the data reduction (Ej2l), the 
revision of the Large catalog (leading to a new version of 
v2.0) and the construction of the new Deep catalog are 
outlined in ^S] and [JH In ^Sl we explain the construction 
of the final Joint catalog which is described in detail in 
ij6] where we also present all the tests and corrections 
involved. A summary is given in ^ 

2. OBSERVATIONS AND DATA REDUCTION 

The central seven pointings of the mosaic for the VLA- 
COSMOS Large project (p ointing numbers: F07 F08, 
Fll, F12, F13, F16, F17; see lSchinnerer et al.l[2007n were 
observed using the VLA A configuration in February and 
March 2006. The observations were executed on 11 days, 
typically lasting for a total of six hours starting at 07hr 
LST (Local Sidereal Time). The exact coordinates of the 
pointing centers are listed in Tab. [1] We used the same 
spectral set-up and set of calibrators as for the VLA- 
COSMOS Large project to allow for an easy combination 
of the data from the two projects: the VLA standard L- 
band continuum frequencies of 1.3649 and 1.4351 GHz 
in multi-channel continuum mode, the quasar 0521-1-166 
(3C 138) for fiux and bandpass calibration at the be- 
ginning of each day, and the quasar 1024-008 as phase 
calibrator. In order to obtain a good sampling of the 
uv-plane, we cycled three times through the pointings 
each day with a total integration time per pointing of 
~45 min. The final resulting integration time per point- 
ing was '-^8.25hr. In addition, we changed the order of 
the observations of the individual pointings each day to 
sample the uv plane uniformly in all pointings. 

During the data reduction process, we excluded all 
EVLA antennas^'' that were present in the array during 
the observations. We followed the dat a reduction pro- 
cedur e adopted for the Large project (jSchinnerer et aLl 
|2007| ) which us es the Astrono mical Imaging Processing 
System (AIPS; iGreisenI 120031 ). uv data points afi'ected 
by RFI (radio frequency interference) were again flagged 
by hand using the task TVFLAG and by excluding ampli- 
tudes greater than 0.55 Jy. The fiux density of the phase 
calibrator 1024-008 was found to be close to its previous 
value in the season of 2004/5. After final calibration, 
the Deep data were combined with the uv data of the 
corresponding pointings from the Large project. 

The final imaging was performed using the task IMAGE 
with the same settings as for the Large project, i.e. 
the same tiling of the individual pointings and the as- 
sociated CLEAN boxes were used, the flux cut-off was 
kept at 45/xJy, and the restoring beam was again set 
to 1.5" X 1.4". The robust parameter was changed from 
to 0.25 to obtain a dirty beam as close as possible to the 

^^ A total of 3 EVLA antennas were included in the array at 
the time of the observations. 



flnal CLEAN beam. As for the Large project, each polar- 
ization for each intermediate frequency (IF) was CLEANed 
separately and combined afterwards in the image plane. 
The final seven new fields were combined with the exist- 
ing 16 fields from the Large project to obtain a mosaic 
covering the fuU 2 deg^ of the COSMOS field. The aver- 
age rms achieved in the central 30' was 8 /iJy/beam at the 
resolution of the Large project. In order to mitigate the 
effect of bandwidth smearing on the derived peak flux 
densities during source extraction (see §6.ip . the final 
mosaic was convolved to a lower resolution of 2.5" (Fig. 
[1]) leading to an increase of the rms to ^12/iJy/beam. 
We tested several convolution kernels by comparing the 
source fluxes and the number of sources extracted for 
two representative sub-images of the mosaic. The res- 
olution of 2.5" provided the best compromise between 
minimizing the flux losses and maximizing the number 
of extractablc sources. We verified that the noise distri- 
bution in the Deep image follows a Gaussian distribution 
as expected (Fig. [2]). This Deep mosaic was used for 
source detection and extraction (see ^. Fig. [3] shows 
the fractional area covered as a function of rms for the 
Deep and Large mosaics. 

3. THE REVISED LARGE PROJECT CATALOG 

The original source cat alog for the VLA-COS MOS 
Large Project presented by iSchinnerer et al.l ()2007D was 
created by the following procedure: 1) the AIPS task 
SAD extracted candidate components down to a peak flux 
density limit of 3cr, 2) after fits to the peaks of these 
candidates using MAXFIT only those with an measured 
(rather than the fitted Gaussian) peak flux density of 
>4.5tT were kept, 3) components potentially missed by 
SAD were recovered with JMFIT which was used to fit 
a Gaussian to all peaks above our detection threshold 
that were not included in the SAD component list, and 
finally 4) components not in the SAD list but with non- 
zero major axes found in step 3 were added. This com- 
ponent catalog was transformed into a source catalog, 
after multi-component sources, e.g. large radio galaxies, 
were identified by visual inspection. Since in the case of 
faint sources, the position and peak values of a Gaus- 
sian parametric fit might not be a good representation of 
the real values, the position and value of the peak found 
by MAXFIT replaced the results from the Gaussian fit. 
Furthermore, all sources classified as unresolved (for our 
methodology see Section 6.2 in Schinnercr et al. 200^) 
had their integrated fiux density set equal to the peak 
flux density. The flnal VLA-COSMOS Large project cat- 
alog (vl.O) has a total of 3601 entries. 

Subsequent use of the Large catalog in conjunction 
with the COSMOS optical and (near- to mid-)IR cat- 
alogs and images showed that the fraction of spurious 
sources increases significantly below 5cr. Spurious sources 
in our definition are sources that lie in a noisier environ- 
ment than expected for a > 4.5ct detection, i.e. three or 
more 3a peaks are present in a 10"xlO" box centered on 
the source. The misidentification is likely due to a mis- 
match in the derived rms value (in a 17.5" x 17.5" box) 
and the local rms distribution at the exact position of 
the source. This finding is also confirmed by a compar- 
ison of sources present in the Large and Deep catalogs 
(sec ii5.2.1[) . Thus we exclude all sources below 5(j from 
the revised Large catalog. In addition, the correction for 
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bandwidth smearing as derived by iBondi et alj (|2008[ ) 
has been apphed to the peak flux densities and the clas- 
sification into resolved and unresolv ed sources has bee n 
modified accordingly (for details sec lBondi et aljr2008f ). 
The revised version of the VLA-COSMOS Large project 
catalog [hereafter: Large catalog (v2.0)] contains a new 
total of 2417 sources. The revised VLA-COSMOS Large 
project catalog (v2.0) has been publicly released to the 
COSMOS archive at IPAC/IRSA". 

4. IDENTIFICATION OF DEEP PROJECT 
SOURCES 

As the rms is less uniform in the Deep mosaic than 
in the Large mosaic and shows a steep increase towards 
the edges, applying a simple flux density cut for source 
detection is not possible. In order to use the AIPS 
source/component finding task SAD, a S/N map was cre- 
ated. First, a sensitivity map (Fig. |4]) was constructed 
using the AIPS task RMSD with a box size of 105" x 105" 
to estimate the rms and a sampling step size of 2.45" in 
both RA and DEC. To obtain the final input S/N map 
for SAD, the Deep image was divided by the sensitivity 
map and blanked at the edges where the rms values in 
the sensitivity map exceeded 34^Jy/beam (the largest 
rms value present in the Large mosaic). This last step 
allowed us to search for sources in the corners/edges of 
the mosaic, whereas the area corresponding to the COS- 
MOS field was searched for the con struction of the Large 
catalog fsee lSchinnerer et al.ll2007l ). 

SAD was set to search for sources in six iterations with 
cut-off levels for the peak of 10, 8, 7, 6, 5, and 4(7. The de- 
rived values in units of S/N were converted to mJy/beam 
using the corresponding rms values from the sensitivity 
map. As in the case of the Large project sources with 
complex structures - e.g. radio galaxies - were fitted 
by multiple Gaussian components, such that the derived 
catalog contains components rather than real sources. In 
order to identify these multi-component sources, we com- 
pared by eye the location of the multi-component sources 
from the Large catalog (v2.0) with the component list 
obtained for the Deep image to identify the components 
belonging to the same source. During the comparison 
with the Large catalog (v2.0), it was noticed that sev- 
eral slightly extended sources were fitted by two separate 
components with a typical separation of ~1" (hereafter: 
'twin' sources). In order to determine the nature (simple 
extended Gaussian vs. truly non-Gaussian geometry) of 
these 'twin' sources, a single Gauss component was fitted 
to these sources using JMFIT. If the integrated flux den- 
sity of this single Gaussian fit was consistent within the 
formal errors with the combined integrated flux densities 
of the two components derived by SAD, we replaced the 
'twin' sources by a single object. 

5. CONSTRUCTION OF THE JOINT SOURCE LIST 

Based on the size d i stribution of faint radio 
sources (e.g. iMuxlow et all l2005t iFqmalont et al.l 120061 : 



lOwen fcMorrisonI 120081 : IBondi et al.ll2008f ) a significant 
fraction of sourc es should be resolved at our fiux limit 
and resolution. lOwen fc MorrisonI ()2008f ) showed that 
source extraction at different resolutions is ideal to max- 

^^ |http : //ir sa ■ ipac ■ caltech ■ edu/data/COSMOS/t ables/vla/ | 



imizc the number of sources found. Therefore, the com- 
bination of the VLA-COSMOS Large and Deep catalogs 
is ideal to find a maximum number of sources for inclu- 
sion in a joint catalog. In the following we describe the 
steps taken to identify such radio sources in the VLA- 
COSMOS mosaics. 

In order to construct a list of radio sources from the 
two separate catalogs, we used the following process that 
is described in detail below: First, sources present in 
both catalogs were included ( §5.2p . then selected sources 
detected only in the Large Project mosaic were added 
( §5.3.ip . as were - finally - a number of sources present 
only in the Deep catalog ( ^5.3.2p . In order to minimize 
the inclusion of spurious sources at the low significance 
end, the local rms was re-examined at 1.5" and 2.5" res- 
olution. All sources having a final peak S/N > 5 at at 
least one resolution were included. 

An overview of the number of sources selected from the 
Large (v2.0) and Deep catalog for inclusion in the Joint 
catalog is given in Tab. [2l together with a listing of some 
key properties. 

5.1. Estimate of the Local Background Noise 

Based on the experience with the Large catalog (vl.O) 
it is important to obtain an accurate estimate of the lo- 
cal rms to properly measure the S/N of individual radio 
sources. In order to find the optimal box size for the 
local rms estimate, we compared the value in the rms 
map used for source detection (see 311) to that obtained 
if a smaller box is used. We chose a box size of 50 pixel 
(=17.5") to obtain a good estimate of the rms around 
compact sources. Fig. [5] shows the different rms values 
obtained with a 105" and 17.5" box for a number of low 
S/N sources from the deep catalog and highlights the im- 
portance of deriving an accurate local rms estimate. The 
choice was based on the finding that a box size of 17.5" 
used for the Large catalog still resulted in a number of 
spurious sources at the low S/N end and the requirement 
that at 2.5" resolution the number of pixels outside the 
source should still be sufficients^. 

We constructed new rms maps at 1.5" and 2.5" resolu- 
tion using a box size of 50 pixels with the task RMSD and 
subsequently used them to estimate the local rms. We 
refer to this revised S/N as the 'local' S/N (S'/7Vi7.5"). 
All sources with a S'/A^i7.5" > 5 will be included in the 
joint source list. 

5.2. Sources Present in Both Catalogs 

For many sources in the Deep catalog it was pos- 
sible to identify counterparts in the Large catalog 

-^^ At an image resolution of 9b (FWHM), a 17.5" rms box holds 

A^bcam = (17-5")7nbca„ ^Vh = (17-5")>p9| X rjh (1) 

independent, hexagon ally packed b eams. Here rji^^O.d is the circle 
packing density (e.g.. rWellslll981 | ) and fcp = 1.13 for an aperture 
with a Gaussian power pattern l lKclll96"^ : IOt oshi 1981). Hence, as 
-'Vbcam(2.5") ~ 40, the 5 a rms estimate has an uncertainty of about 
5 (t/\/40 Ri 0.8 a, which implies that the local S/N could be 4.2 a or 
less in some regions of the mosaic. To get an upper bound on how 
many of our sources that were only detected with S/Nijr^n > 5 at a 
resolution of 2.5" might be affected by this, we checked which also 
had S/N > 5 in 105" rms boxes and found that 61% also satisfied 
the latter criterion. These were assigned a flag 'det' = 1 in the Joint 
cata log, while the remaining sources are given 'det' = 2, see Section 
[0]and Tab. |3] 
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(jSchinnerer et al.]|2007D down to 5a or even 4.5cr. The 
successfully matched sources include all but one of the 
objects already known from the Large Project to con- 
sist of multiple Gaussian components and the majority 
(^60%) of the single component sources in the Large 
catalog (vl.O). 

5.2.1. Identification of Single-Component Sources 

To identify radio sources present in both the Deep and 
Large catalog, we used the Large catalog down to 4.5(t. 
We matched the two catalogs with a search radius of 1.5" 
to account for offsets between the peak positions for ex- 
tended objects. The spatial characteristics of the noise in 
the Deep and Large image differ due the additional cov- 
erage of the central seven pointings by the Deep project 
and due to the differences in the final resolution (2.5" in 
the Deep mosaic vs. ^1.5" in the Large mosaic; see ij2|). 

As mentioned in Section [3] the occurrence of spurious 
sources in the Large catalog is significantly reduced by 
omitting sources in the range S/Ni^arge G [4.5, 5[. The 
added depth gained with the Deep project observations 
in the central 40' (Fig. |6|) should increase the reliabil- 
ity of sources with 5/A^Largc < 5 in the Large Project 
catalog (vl.O) within this area. Fig. [6] shows the cor- 
responding detection thresholds of the two images. The 
peak flux densities of sources in the Large catalog (vl.O) 
are plotted as a function of their radial distance from 
the field center, their lower envelope roughly tracing the 
average S/N detection threshold of 4.5cr (used as cut for 
the Large catalog (vl.O)). The radial evolution of the 
mean rms (rms) was calculated in concentric rings with 
a width of 2' using the respective sensitivity maps for 
both images (Fig . Ill for the Deep mosaic and Fig. 12 in 
iSchinnerer et al.r(|2007[ ) for the Large one). The 'bump' 
in the measured mean noise level at a radius r « 45' is 
due to the sensitivity maps' rectangular geometry. In or- 
der to obtain an estimate of the minimal average noise 
level at this distance from the field centre we fitted a 
cubic spline to the measured rms in this region. The re- 
sulting curve for the Large mosaic, when scaled to 4.5cr, 
defines a lower envelope to the measurements (ignoring 
detections in sites with locally lower rms noise). The 4cr 
limit of the Deep catalog corresponds to the 4.5cr limit of 
the Large catalog for the inner radii with r < 15' while 
it increases to the 5ct limit of the Large catalog at larger 
radial distances. 

The fraction of sources in the VLA-COSMOS Large 
catalog (vl.O) which are at a given distance from the 
field centre and have S'peak < 4(TDeep is shown in Fig. |B]d. 
The increase of this fraction with distance from the field 
center will affect the radial dependence of the number of 
successful matches between the two catalogs. We veri- 
fied that the success rate of the catalog matching shows 
no distance dependence (Fig. [7]). While the fraction 
of sources in the Large catalog (vl.O) for which a coun- 
terpart could be identified decreases at larger distances 
(solid black line), the effective matching success rate de- 
fined as ratio of the measured matching success rate and 
the fraction of sources with peak flux densities below Aa 
in the Deep mosaic stays basically constant (thick grey 
line) within a range of 60 and 80%. 

Thus we conclude that the direct positional correlation 
of catalog entries has not caused sources to be system- 
atically missed as a function of (radial) field distance in 



the VLA-COSMOS Deep mosaic. 

Before final inclusion in the source list we checked for 
each of the successfully matched sources that it met our 
selection criterion of peak S/Nn,^// > 5 at either 1.5" or 
2.5" resolution in the deep data. If this requirement was 
not satisfied the corresponding source was not considered 
any further. 

5.2.2. Identification of New Multi- Component Sources 

As already stated above, all but one of the multi- 
component sources listed in the Large catalog (v2.0) were 
also present in the Deep catalog. However, during the 
merging of the two catalogs 43 new multi-component 
sources were identified. Image cutouts of all new multi- 
component sources are shown in Fig. [51 

During the cross-correlation of the Large and Deep 
catalog several sources were recognized to constitute 
sources made up of several (Deep) rather than just a 
single (Large) Gaussian component during visual inspec- 
tion of all sources. In addition to this, several Deep 
sources found in the neighborhood of multi-component 
objects (identified in the Large catalog) were subse- 
quently assigned to these sources as additional com- 
ponents of their extended emission. In one instance 
two of the former multiple component sources were 
joined to form a new, larger multiple component source 
(a bipolar jet with nucleus assigned the new Joint 
catalog ID COSMOSVLADP_J095758.044-015825.2: - 
the three merged components have the following IDs 
in the VLA-COSMOS Large Project catalogs: COS- 
MOSVLA.J095755.84+015804.2, J095758.04-1-015825.2 
& J095800. 79+015857.1). Besides these newly identified 
or augmented multi-component sources there was also a 
small number of multi-component sources that had ei- 
ther (i) not been previously listed in the Large catalog 
or (ii) are situated outside the area searched during the 
construction of the Large catalog. 

5.3. Sources present in only one catalog 

Apart from the sources that are present in both cata- 
logs, objects present only in one of the two catalogs have 
been added to the list as well. The exact procedures 
adopted are described in the following. 

5.3.1. Sources identified m the Large image only 

All sources in the Large catalog (v2.0) with a 
S/Ni^^^gc > 5 but without a match in the Deep cata- 
log are in principle valid candidates for inclusion in the 
joint source list, provided they are not flagged as detec- 
tions potentially due to side lobes of strong radio sources. 
Since SAD tries to fit Gaussians to all peaks above a cer- 
tain limit, the fact that certain sources were not found by 
SAD during the construction of the Deep catalog suggests 
that some of these objects might be mere noise peaks 
or have very unusual morphologies. Similarly, we expect 
that a fraction of the sources with S/N < 4.5 might be 
real. 

Therefore, we measured for all sources in the Large 
catalog (vl.O) without a match in the Deep catalog the 
peak flux density in the Deep image at 1.5" and 2.5" 
resolution at the according position from the Large cat- 
alog following the steps described in ij6.2.1l and subse- 
quently discarded 959 candidates with S/Nn,^" < 5 at 
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both resolutions while it was possible to keep 393 of the 
unmatched sources with S'/iVLargo > 4.5 from the Large 
catalog (vl.O). 

5.3.2. Sources identified in the Deep image only 

1155 objects - predominantly in the range of 4< 
S/Nueep <8 ~ were detected in the Deep mosaic which 
did not have a counterpart in the original Large catalog 
(vl.O). In order to assess the significance of these sources 
which arc potentially new detections due to the different 
source extraction area for the Deep catalog, the increased 
sensitivity in the central square degree and the fact that 
the larger beam is more sensitive to slightly extended 
sources (as SAD works with peak fluxes), we adopted the 
approach outlined below. While the S/Nuccp = 4 cut 
corresponds at least to the detection threshold used for 
the original Large catalog, our experience with this cat- 
alog showed that the number of spurious sources rises 
significantly below S/N < 5 due to small mismatches be- 
tween the real local rms and the estimated value in the 
constructed rms map. 

Thus, the local S/N was checked at 2.5" resolution as 
well as at 1.5" resolution. All sources with S/Nn,5» > 
5 were included in the joint list. In total, 178 of 283 
unmatched Deep detections with S'/iVDeep ^ 5 were 
taken over in the Joint catalog. Of 872 sources with 
4 < S/N-Dccp < 5 some 136 met our selection criteria. 

6. JOINT CATALOG 

In the following we describe how the flux densities and 
all other measurements provided in the catalog were de- 
rived and which corrections have been applied to the tab- 
ulated values. 



6.1. Correction for Bandwidth Smearing 

Bandwidth smearing (BWS) can occur in radio synthe- 
sis imaging when a finite frequency range is observed. It 
causes a radial smearing that becomes more severe with 
increasing distance from the phase center (or center of 
a pointing), as the phase calibration is mathematically 
speaking only correct for a given frequency at the phase 
center. Thus the effect is similar to chromatic aberration 
in optical imaging. A detailed explanation and discus- 
sion of the effe ct of bandw i dth sm eari ng for the VLA is 
;iven e.g., by iThompsoiil (J1999D and IBridle fc SchwabI 
19991) . 

For the VLA-COSMOS project with an observed band- 
width of 3.125 MHz for a single channel, a (final) radial 
smearing of 2.25" is expected at a radial distance of 15' 
from the pointing center, i.e. close to the cut-off radius of 
the individual pointings used when creating the mosai c 
(for equations and details, see IBridle fc SchwabI 1 1999( ). 
This bandwidth smearing effect also causes a decrease of 
the peak flux density as the emission is now distributed 
over a larger area. The integrated fiux density, however, 
will not be affected. As sources in the final mosaic can 
be covered by up to 7 individual pointings, the impact of 
bandwidth smearing strongly varies at each location in 
the mosaic as a given source is separated by a different 
amount from the centers of the individual pointings. 

Thus, the effect of bandwidth smearing depends on 
the resolution of the image. In correcting the measure- 
ments presented in the Large catalog (v2.0) for the BWS 



effect iBondi et al.l (|2008D adopted the approach of boost- 
ing the measured values of the peak flux density Speak 
by different amounts depending on the distance from the 
center of the VLA-COSMOS field. The boost factors 
were derived from tests comparing the peak fiux den- 
sities of sources at the centers in individual pointings 
and the derived flux densities in the final Large mosaic. 
This approach c an be used becau s e the sensitivity map 
(cf. Fig. 12 in ISchinnerer et al.l (120071) ') for the VLA- 
COSMOS Large project is fairly uniform over significant 
areas, i.e. all locations suffer similarly from the combined 
BWS effect of the contributing pointings. However, the 
additional data for the central seven pointings for the 
Deep Project changed the uniformity of the sensitivity 
distribution over the COSMOS field. W e thus exam- 
ine both the method of IBondi et al.l (J2008D as well as an 
alternative based on a modeled sensitivity map for the 
VLA-COSMOS Deep project. 

The model sensitivity distribution for an individual 
pointing was constructed using the beam pattern for a 
VLA antenna at 1.465 GHz, as specified in the help page 
for the MPS task PBCDR: 



1-1.343x10"^ j:'2+6. 579x10"^ JT^'-l. 186x10"^° JT'^ , 

.(2) 
where ^ is the product of the distance from the point- 
ing center in arcminutes with the observing frequency 
in GHz. The full model sensitivity map (Fig. [TT|) was 
then assembled by adding the individual beam patterns 
which were weighted according to the observation time 
dedicated to each pointing (i.e., the central 7 pointings 
were weighted more strongly than the others by a factor 
of 1.4). 

To estimate the magnitude of the required BWS cor- 
rection we c losely follow the procedure described by 
IBondi et al.l (|2008f ). In brief, a primary beam correction 
was applied to the images of the 23 individual point- 
ings and they were convolved to the same resolution of 
FWHM 2.5" as the Deep mosaic. Then we measured 
the peak flux densities of selected sources in these in- 
dividual images and compared their values with those 
obtained from the same measurement carried out at the 
corresponding position in the mosaic. 

At this stage, it is mandatory to use sources with a 
peak flux density that has not been signiflcantly dimin- 
ished by the BWS effect. The dimensionless parameter 



/3 



Ai. 



fsee lBridle fc SchwabI ([T999h and the VLA 



1^0 fpWHM 

observational status summary^'^) makes it possible to in- 
fer the amount by which peak flux densities are reduced. 
In the case of the VLA-COSMOS Large image with an 
adopted beam width of ^fwhm ^ 1-5" the peak flux den- 
sities of sources within 5' of the pointing centers are ex- 
pected to be reduced by less than 5%. At a resolution of 
^FWHM ~ 2.5", as in the case of the Deep image, the same 
small decrease still exists as far as 8' from the pointing 
center. Thus we were able to use a contiguous area cov- 
ering a significant fraction of the entire VLA-COSMOS 
mosaic (see area indicated by dashed white circles in Fig. 
[TT|) . A correspondingly larger number of sources could 
thus be used to check the validity of the adopted method 
for BWS correction. In the following we will quantify the 

^•^ http : //www . via . nrao . edu/astro/guides/vlas/current/nodelS . html 
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magnitude of the BWS effeet using only sources within 
5' of the center of either of the 23 VLA-COSMOS point- 
ings ~ this ensures that peak flux densities measured on 
individual pointings should differ by less than 2% from 
the nominal value - and we will use the larger number 
of sources out to 8' to eheek the quality of the approach. 
Since the aim is to correct for the BWS effect introduced 
by overlapping pointings, strictly speaking no (valid) cor- 
rections can be derived for sources located in the edges 
(i.e. those covered by only single pointing, see Fig. [T^ . 
Fig. [T2] shows the comparison of the two different 
methods to correct for the BWS effect. It should be 
noted that the maximum decrease observed is only 10%. 
The top row displays the uncorrected values of the quan- 
tity y = Speak, mosaic/S'poak, pointing, 1-6. the ratio of the 

peak flux densities measured in the mosaic and in an in- 
dividual pointing, as a function of normalized sensitivity 
(left column) and radial distance from the field center 
(right column). The median ratio (^) of peak flux den- 
sities in a specific bin is marked by an open circle. The 
associated error bars span the interquartile range of -S^ 
within a bin in distance or model sensitivity. Median and 
errors are reported in red when only sources within 5' of 
the center of one of the 23 pointings are used (i.e. only 
the red points) and in black when all sources within 8' 
are considered (i.e. the red and grey points). 

It is obvious from the upper row of Fig. [12] that the 
effect of bandwidth smearing disappears^^ (i.e. the value 
of the ratio ,5^ approaches unity) towards the edge of 
the field where there are no overlapping pointings or, 
conversely, at lower sensitivity values which stem from 
regions in the field covered by only one pointing. The 
best-fitting linear trends define which BWS correction 
should be applied to sources lying at a given distance or 
within a region of a specific sensitivity value, depending 
on the method considered. In the middle row the ratios 
-^sensitivity correction are Corrected according to the second 
method using the model sensitivity map (Fig. [TT|) . i.e. 
the the slope of the line as a function of sensitivity (up- 
per left panel). Thus the distribution of ,5^ is necessarily 
flat when plotted as a function of model sensitivity. This 
flatness is also preserved when the corrected peak flux 
densities are plotted against distance from the center of 
the field. The BWS correction (third row) based on the 
comparison of the peak flux density ratios as a function 
of distan ce from th e center of the VLA-COSMOS field 
following lBondi eta l. (2008) also fares well in straighten- 
ing the distribution of corrected ratios ^distance correction 
with respect to both distance (again, a requirement) and 
model sensitivity. However, a weak correlation of peak 
fiux density ratios with sensitivity is still present after 
the application of this method. Therefore we conclude 
that the sensitivity-based approach is slightly better for 
BWS correction and will use it to compute the corrected 
peak flux density, SV,eak, corr. , provided in the Joint cata- 
log (see also Tab. [3|. 

6.2. Catalog Entries 
6.2.1. Single- component sources 

95% of the sources in the VLA-COSMOS Joint cat- 
alog are well fit by a single Gaussian component. The 

^'^ Note that at these radial distance the BWS effect can no 
longer be determined due to the lack of overlapping fields 



following five steps were required to obtain the position, 
peak and integrated flux density, as well as the size for 
each source. All measurements described in the following 
were carried out on the Deep Mosaic. In particular, this 
also applies to those sources which were not part of the 
initial Deep catalog (see §5.3.1[) . 

1. The (parametric) integrated fiux Stotai was derived 
using the task SAD. In the case of the 'twin' sources, 
JMFIT was used instead (for details see g]). At the 
same time the convolved shape parameters were 
determined. 

2. The peak flux density was obtained using the task 
MAXFIT^^. The square search box for MAXFIT was 
centered on the position of the Gaussian fit from 
the previous step. We used an initial box size of 
3 pixels that was increased in steps of 2 pixels if 
MAXFIT did not converge. The measured peak flux 
density S'pcak was corrected for the effects of band- 
width smearing (see §6.1|) resulting in the corrected 
measured peak flux density S'pcak, corr.- 

3. The source position was set to the position of the 
MAXFIT peak. 

4. As a significant fraction of sources is resolved at 
the faint fiux levels probed here, it is necessary to 
classify the sources into resolved and unresolved. 
We fo llowed the approach used bv lSchinnerer et ahl 
()2007[ ) anddondi et al. (2008). The method relies 
on the fact that the ratio between integrated and 
peak fiux density is a measure of the spatial extent 
of a radio source. A detailed description of the 
method is given in Appendix |A] together with a 
description of a suite of Monte Carlo simulations 
which we use to calibrate our classification scheme. 
Fig. [13] is the diagnostic diagram used to identify 
the resolved sources based on our simulations. The 
line defining the lower boundary to the locus of 
resolved sources in Fig. [13] is given by the following 
equation^^: 



Stotal/ Speak, corr. — U.OO 



-iy(s/iv)i-*^ 



(3) 



A discussion on possible systematic biases inherent 
to this method is presented in the Appendix. We 
performed an independent check of this classifica- 
tion me thod by adopting th e JM FIT approach used 
by, e.g ., lMriler et alJ () 2008i) and lOwen fc MorrisonI 
()2008f ). As JMFIT also derives upper and lower lim- 
its on the deconvolved size of the major and minor 

^^ Note that MAXFIT derives the maximum by fitting a quadratic 
function to a 3x3 pixel map. Given that our Large (Deep) CLEAN 
beam is well sampled writh at least 4 (7) pixels per axis the map 
value and the MAXFIT value agree within ~0.2%. 

^^ Note that in previously published catalogs (see references 
above) using this approach, the line described by eq. ^ follows 
from logarithmically mirroring a lower envelope to the points bel ow 
the line S'total/'S'pcak.corr. =1 abovc this very line. In Fig. 1131 we 
show that the negative mirror image of eq. Q is a very conserva- 
tive lower envelope of the kind just described. As a consequence, 
the classification into resolved/unresolved sources adopted in this 
paper based on the simulations in the Appendix, leads to a signif- 
icantly larger fraction of sources being classified as 'unresolved' in 
the Joint catalog. 
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axis of a source, these limits can also be used to 
divide sources into resolved and unresolved. If the 
upper limit of the major axis was equal to zero, a 
source is considered unresolved. Comparison be- 
tween both methods shows that <1% of our re- 
solved sources would be considered unresolved us- 
ing the JMFIT criterion while 66% of the unresolved 
sources would be resolved. Thus the total fraction 
of resolved sources would i ncrease to ~ 7 0% si milar 
to the fractio ns found W iMiller et all ()2008[ ) and 
lOwen fc Morr ison (200^. Using equation (jS]) we 
have found 405 resolved sources for which the inte- 
grated flux is given by the total flux of the Gaus- 
sian sad/ JMFIT fit, and 2329 unresolved sources for 
which the integrated flux is a posteriori set equal 
to the corrected peak flux density. These num- 
bers correspond to 14% and 81% of the sources in 
the Joint catalog, respectively, while the remaining 
5% represent radio sources with multiple Gaussian 
components (see following section). 

5. The deconvolved source size parameters (major 
axis OM,dec, minor axis 0m, dec and position angle 
PAdoc) were calculated for the resolved sources ac- 
cording to the algorithms of the MPS task JMFIT. 
For unresolved sources all three are set to zero. 

6.2.2. Multi-Component Sources 

For sources consisting of multiple Gaussian compo- 
nents the integrated flux density was manually measured 
on the VLA-COSMOS Deep image with the AIPS task 
TVSTAT in order to encompass all the emission from these 
irregularly shaped sources. As the individual compo- 
nents all have flux density peaks of their own, no peak 
value is specified for multi-component sources in Tab. [3l 
Positions were adopted from the VLA-COSMOS Large 
project catalogs whenever possible. In all other cases the 
position was chosen on an individual basis based on the 
radio morphology of the sources. In practice this usu- 
ally corresponds to a luminosity weighted mean of the 
positions of the subcomponents (as already done for the 
multi-co mponent sources in th e Large Project catalog, 
see §6 of lSchinnerer et ani2007| ). 

Measurements of the major and minor axes along the 
direction of maximal extension of the source and at right 
angles to the latter, respectively, are provided for multi- 
component sources. However, no estimate of the position 
angle is provided as it is usually impossible to define for 
the complex geometry of these sources. Note that while 
the coordinates were taken over, fluxes and dimensions of 
multi-component objects adopted from the Large catalog 
(v2.0) were re-measured in the Deep project mosaic at 
2.5" resolution. 

6.3. Errors on the Catalog Entries 

Below we summarize how the errors on the source prop- 
erties listed in the VLA-COSMOS Joint catalog were 
computed. In the case of the single-component ob- 
jects analytic formulae from the literature are adopted. 
For the multi-component objects which have a manually 
measured integrated flux density S'totai, the accuracy is 
about 10% based on the comparison of repeated flux mea- 
surements on the same objects by different people. For 
multi-component objects no error estimates are provided 



for the source shape parameters as these are meant as a 
rough measure of the angular size of the source. For the 
same reason we also do not provide values for the er- 
ror on the major and minor axes of the Gaussian single 
component sources. 

The error on the peak flux density of single component 
sources is defined as the local rms noise at the position 
of the source. For unresolved sources the integrated flux 
density and peak flux density are identical and hence the 
error on the integrated flux density also corresponds to 
the local rms noise. 

In calculating the error crstotai ^r resolved sources 
in the VL A-COSMOS Jo i nt ca talog we fo llow the ap- 
proach of iHopkins et al.l (I2003D (see also iBondi et al.l 
l2003t iSchinnerer et al.l 120041 ) based on the assumption 
that the relative uncertainty crg^^j^j/S'totai in the inte- 
grated flux density is due to uncertainties /idata in the 
data and uncertainties /ifit in the Gaussian fit: 



O'Stotal 

•Stotal 



Y Mdata Y 
V '5'total / 



V "Stotal / 



(4) 



The two factors beneath the squ are r oot in th e equa tion 
are fsee [Windhorst et al.l (I1984D andlCondonI (I1997D . as 
well as the explanations in ISchinnerer et al.l (120041 )): 



(where S/N 



Mdata 
'S'total 

= -Spcak/rms) and 




0.012 



Mfit 



s, 



total 




(5) 



— + ^ . (6) 



Here 9b and 9b are the major and minor axis of the beam 
(in the present case of a circular beam 9b = 9^ = 2.5") 
and, analogously, 9m and 9^^ the major and minor axis 
of the measured (i.e. convolved) flux distribution. The 
S/N values of the fit - ps, pM, and pm - arc parameter- 
dependent: 



Px 



7Mfm 



9m I \ 9r. I \ N 



(7) 



^9B9b 

with a = j3 = 1.5 for ps, a ~ 2.5 and /3 — 0.53 for pM, 
and a = 0.5 and /3 = 2.5 for p,u- 
The errors on the position of single-component sources 



are 




4 ln(2)p| 



4 ln(2)p2 



sin(PA)" 



where €a ~ ^s — ^b/IO ~ 0.25" are the positional 
calibration errors on right ascension and declination and 
PA is the measured position angle of the Gaussian fit. 

6.4. Description of the VLA-COSMOS Joint Catalog 
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The information on the procedures and equations used 
to calculated each entry in the VLA-COSMOS Joint 
catalog has been presented above. For each source we 
list the source name, the ID of the source in the VLA- 
COSMOS Large Project catalogs (where available), as 
well as the derived source properties and the associated 
errors. Furthermore, the source properties and flags 
associated with each object are explained. All 2865 
radio sources in the VLA-COSMOS Joint catalog are 
listed by increasing right ascension in Tab. [3] with the 
following columns: 

Colunin(l): Source name 

Column(2): Source name in VLA-COSMOS Large 
Project catalog (set to J999999.99+999999.9 if not 
present in the VLA-COSMOS Large catalog) 



Column(3) 

Column(4) 

Column(5) 
format 



Right ascension (J2000.0) in degrees 
Declination (J2000.0) in degrees 
Right ascension (J2000.0) in hexagesimal 



Column(6): Declination (J2000.0) in hexagesimal format 

Colunin(7): rms uncertainty in right ascension in arc- 
seconds 

Colunin(8): rms uncertainty in declination in arcseconds 

Colunin(9): Peak flux density and its rms uncertainty 
in mjy/beam 

Column(lO): Peak flux density corrected for bandwidth 
smearing in mJy/beam 

Column(ll): Integrated flux density and its rms uncer- 
tainty in mjy 

Column(12): rms measured in the RMSD sensitivity map 
in mJy/beam 

Colunin(13): Deconvolved major axis size 0m, doc in 
arcseconds 

Column(14): Deconvolved minor axis size 0m, doc in 
arcseconds 

Column(15): Deconvolved position angle of source 
PAdoc (counterclockwise from North) in degrees 

Column(16): Flag for resolved and unresolved sources 
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-2 - unresolved only in VLA-COSMOS Deep image 
-1 - resolved only in VLA-COSMOS Large image 

- unresolved in both VLA-COSMOS Large & 

Deep images 

1 - resolved in both VLA-COSMOS Large & Deep 

images 

2 - resolved only in VLA-COSMOS Deep image 

Column(17): Flag for distinction of multi-component 
and single-component sources 

'^^ Note that the flag value of '-2' is only assigned to sources 
that were not in the VLA-COSMOS Large Project catalog. The 
flag '-!', on the other hand, is used for sources that were classified 
as resolved in the Large Project catalog but are listed as unresolved 
in the Joint catalog. 



- single component source 

1 - multi-component source identified in VLA-COS- 

MOS Large image 

2 - multi-component source identified in VLA-COS- 

MOS Deep image 

Column(18): Flag for catalog membership 
-1 - only detected in the VLA-COSMOS Large image 

- detected in both the VLA-COSMOS Large & 

Deep image 

1 - only detected in the VLA-COSMOS Deep image 

Column(19): Flag specifying at which resolu- 
tion the source was detected with S/Nn,^n > 5 
-1 - detected with S/Nn 5/' > 5 only at a resolution 
of 1.5" 

- detected with S/Nn,z» > 5 at both 1.5" and 2.5" 

resolution 

1 - detected with S/N > 5 only at a resolution of 

2.5", but in both the large and 

small scale (105" and 17.5" box size, respectively) 

rms map 

2 - detected only with S/Nn,5>> > 5 (but not in the 

large scale rms map) at a 
resolution of 2.5" 



6.5. Comparison between Large (v2.0) and Joint 
Catalog 

We compared the source properties in the flux-limited 
Large catalog (v2.0) and in the Joint catalog for those 
2256 catalog entries which are common to both catalogs. 
The different resolution in the VLA-COSMOS Deep and 
Large mosaics leads to a different sensitivity to extended 
emission. Hence, source shape parameters cannot be ex- 
pected to show exact agreement between the two cata- 
logs. Therefore the comparison is limited to the peak 
and integrated flux densities. The results of the compar- 
ison are illustrated in Fig. [15] and Fig. [16] for the peak 
and integrated flux density, respectively. In particular, 
it is instructive to compare the peak values of sources 
that are classifled as unresolved in both catalogs. These 
sources are highlighted in light grey in Fig. [15] and are 
seen to scatter well around the diagonal line of 1:1 cor- 
respondence between the measurements from the Joint 
and Large mosaics. The comparison of the integrated 
flux densities (Fig. [T6| shows that there is a small popu- 
lation of objects (light grey points) lying above the diago- 
nal (i.e. their integrated flux density is higher in the Joint 
catalog). These objects are sources that were previously 
classified as unresolved but have acquired the status of re- 
solved in the Joint catalog, implying that their integrated 
flux densities have necessarily increased when measured 
in the Deep mosaic. This class of objects accounts for 
nearly all of the outliers in the plot. Multi-component 
sources (dark grey crosses) in general lie also above the 
diagonal bisector, reflecting the enhanced sensitivity to 
extended emission in the Deep image. For most other 
objects the new and previous estimate of integrated flux 
density is in good agreement and hence also with flux 
density e stimates of the NVSS and FIRST surveys as 
shown in iSchinnerer et ahl (|2007[ ). 
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Contin ued analysis of th e VLA - COSMOS catalog pre- 
sented in iSchinnerer et all (|2007[ ) and the completion of 
the VLA-COSMOS Deep project motivated the compi- 
lation of a new radio catalog for the COSMOS field. 
The VLA-COSMOS Joint catalog was generated by com- 
bining the catalogs of the VLA-COSMOS Large Project 
with a newly created source catalog (Deep catalog) from 
the 2.5" resolution Deep mosaic. This catalog is already 
available for download by the public at the COSMOS 
archive at IPAC/IRSA^^. A comparison of the depth and 
areal coverage for a representative sample of deep field 
radio surveys at 1.4 GHz shows that the VLA-COSMOS 
covers the largest area at its depth and angular resolu- 
tion (Fig. [IT)) . Thus it should be well suited to also 
study effects of the Large Scale Structure on the pres- 
ence/absence of radio emission. 

The reduction and analysis of the deeper 20 cm obser- 
vations of the central 7 pointings of the VLA-COSMOS 
projects using the VLA in A configuration have been de- 
scribed in detail (also referred to as VLA-COSMOS Deep 
project). In order to minimize the effect of bandwidth 



smearing the Deep mosaic has a resolution of 2.5" (com- 
pared to 1.5" for the Large mosaic) and it was used to 
create a corresponding source catalog using the task SAD. 

An input list for the new Joint catalog was compiled 
by combining the revised Large catalog (v2.0) and the 
new Deep catalog. The criteria were set such that no 
particular bias against slightly extended radio sources 
was present when selecting the sources. All properties 
of the radio sources listed in the Joint catalog have been 
derived in the 2.5" resolution Deep mosaic. 

The construction of the Joint catalog was motivated 
by the desire to provide a catalog of bona-fidc radio 
sources in the COSMOS field for distinct science appli- 
cations that are interested in the radio properties of cer- 
tain populations of galaxies. On the other hand the re- 
vised Large catalog (v2.0) is flux-limited (in radio), has a 
fairly uniform sensitivity coverage and its c ompleteness 
is well characterized (see iBondi et al.ll2008[ ). thus it is 
well suited for , e.g., studies of the faint r adio population 
(such as, e.g.. lSmolcic"eraII[200l l2009al lbl). 



APPENDIX 
BIAS ON THE ESTIMATION OF INTEGRATED FLUXES AND SOURCE SIZES 

At low signal-to-noise, source fitting algorithms (e.g. JMFIT) are increasingly liable to return false results, an 
effect which is further amplified by the fact that there are multiple free parameters (e.g. total flux, two angular 
size components) which cannot be legitimately fixed to make the fit more robust. This Appendix discusses biases 
in the total flux measurements in our catalog which may arise as a consequence of this. We will not study the 
uncertainties, systematic or not, which might affect the angular source sizes we quote in the catalog, as we consider 
these less important for most users of the Joint catalog. Since the unresolved sources in the catalog have a differently 
determined (MAXFIT) total flux, any biases in the Gaussian source fitting will only be of importance for sources classified 
as resolved. How many sources are affected is determined by the choice of classification scheme (see Section [6. 2. ip . 
Part of this section hence deals with our choice of a criterion that both limits fiux biases in the final sample and also 
performs satisfactorily when it comes to separating resolved from unresolved objects. 

We first des cribe our s imulations. We performed a set of Monte Carlo simulations following the approach used for the 
Large project (jBondi et al. 2008). Circular mock sources with random fiux densities down to SOyuJy and si zes following 
the m easured radio source counts and an angular size distribution {6) ^^ S™ with an exponent m=0.5 (see IBondi et al.l 
120081 for details) were inserted into the Deep mosaic. They are not subject to bandwidth smearing effects. We 
searched at the positions of these 16,000 sources using MAXFIT followed by the application of JMFIT (single component 
fit). We verified that this approach gives the same results as using SAD. Roughly 6,500 sources could be recovered with 
a S/N>5. We also ran simulations with different source size distributions which all returned qualitatively similar 
results. We would like to remind the reader that the final numbers and errors quoted below are sensitive to the adopted 
intrinsic angular size distribution and therefore also to the resolution of the mosaic used. 

As a significant fraction of sources is resolved at the faint flux levels probed by the VLA -COSMOS project, it is 



nece ssary to class i fy the sources into resolved and unresolved. We ado pt the approach used by ISchinnerer et al.' (|2007f ) 
and IBondi et al.l (|2008f ) (for applications to other radio surveys, see iPrandoni et all 120001 : [Bondi et al.l |2003() . The 



method relies on the fact that the ratio between integrated and peak flux density is a measure of the spatial extent 
of a radio source (given by the major and minor axis ^m and 0m) hi comparison to the size of the synthesized beam 
(with major and minor axis 0b and 9^,): 



<S'total/5'poak = (^M Om) / {Ob ^b) 



(Al) 



On the other hand we can estimate the limiting (intrinsic) size ^limit a t a given S/N above which a source could be 
classified as resolved. The threshold ^umit is estimated using cq. (16) of lCondonl (jl997| ): according to this expression 
the error in size (A^src.) is proportional to (S/N)~^ (sec Fig. [THl dashed line) and a point source would thus on average 
be liable to have a convolved size ^obs. = Ob^ A^src- AH sources with intrinsic sizes below the lower (red) line in Fig. 
[TBI i.e. ^limit = \/&ths. ~ ^B ~ t^^ inferred size of a point source subject to the S'/TV-dependent error /S.6src., cannot be 
expected to be resolvable. 

The ratio 5'totai/S'poak of our 'unresolvable' (black; 9 < ^umit) and 'resolvable' sources (grey; 9 > 6'iimit) is plotted 
as a function of the S/N with which the simulated sources are detected in Fig. [121 After calibration with the 
simulated sources, this diagnostic diagram will be used to identify the unresolved sources in the real catalog. By 
necessity, measured values with S'totai/'S'poak < 1 are due to the infiucncc of image noise on the determination of source 
flux density and/or size. In general, the noise can both lead to an artiflcial reduction or increase of the true flux 
density ratio (as peak and integrated flux density are determined independently) implying that not all sources with 
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5'totai/'5'pcak > 1 are genuinely resolved either. As can be seen in Fig. [12] there is a fairly well- defined loeus above 
which only few unresolvable sources occur. We approximate it by a line with functional form S'totai/'S'pcak = « 



_ -VCS/N)"^ 



As our main goal is to minimize total flux density biases, we define a rather conservative separating line: 



S'total/'S'pcak — 0.35 



VCS/N) 



(A2) 

This choice ensures that a minimal number of unresolvable moc k sou rces is classified as resolved (in the following a 
'resolved' source is understood to lie above the line given by eq. (|A2[) in Fig. ^^, especially at low S/N where the 
errors on total flux measurements from Gaussian fits are largest. A curve rising more slowly toward low S/N would 
have raised the number of resolvable sources actually classified as resolved. It would also, however, have increased the 
number of misclassified unresolvable sources. 

In the upper window of Fig. [12] we illustrate which fraction of resolvable (unresolvable) mock sources is classified as 
unresolved (resolved) in a given range of S/N. In total, ^38% of the resolvable and ~3% of the unresolvable sources 
are 'misclassified', leading to an overall success rate of nearly 87%. 

Given the classification of our mock sources into resolved and unresolved following eq. (jA2[) . their total fluxes can 
now be set to S'totai == ■S'.jmf it a nd S'totai = -Speak, respectively (see Section [6.2. ip . and any systematic flux biases at 
low S/N quantified. In Fig. 1201 we show the difference between the input {St-m) and output (Stout) integrated flux 
densities as a function of S/N (left-hand column) and as histograms (right). As the mock sources were randomly 
injected into the mosaic and no attempt has been made to avoid confusion with other sources, the outliers in the 
distributions can be explained due to confusion with real sources. However, the median derived from the distribution 
should be unaffected by these outliers. The median relative offset between flux values in different ranges of S/N is 
listed in Tab. [5j In our lowest S/N bin (5 < S/N < 6) we detect a small bias of '^5% (in the sense that recovered 
fluxes overestimate the true value) for all sources which becomes even more negligible for higher significance sources. 
Although this bias is signficantly larger for resolved sources, its effect on statistical tools, such as source counts, is 
minimal, as these use all sources and the fraction of resolved source is at low S/N is very small. As the scatter in the 
derived properties (at fixed S/N) is much larger, we conclude that a systematic correction of the derived integrated 
flux on a source by source basis is neither straightforward nor indispensable. 
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Fig. 1.— The VLA-COSMOS Deep image at 2.5" resolution. 
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Fig. 2. — Distribution of the noise in the VLA-COSMOS Deep mosaic. Pixel values extracted from a 52' X 52' box centered on the 
COSMOS field center show a Gaussian distribution in agreement with our assumption of Gaussian noise. The fitted Gaussian (dashed line) 
has a rms of 12.09 /xjy/beam (ct). The tail at high flux densities is consistent with the presence of sources in the image which were not 
excluded during the extraction process. 
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Fig. 4. — Sensitivity map of the area covered by the VLA Deep Project, derived using the AIPS task RMSD (see text for details). The rms 
is fairly uniform except for areas around strong radio sources. Lighter shades indicate lower rms noise values. The contours correspond to 
rms levels of 10, 12, 15, 20, 25, 30, 34 and 40^Jy/beain. The search for radio sources was limited to the area enclosed by the (light grey) 
34/iJy/beam contour. For comparison, the dashed box indicates the area used to construct the catalog of the Large Project. 
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Fig. 5. — Comparison of the noise estimate in the rms map used to generate the S/N map for source detection in the 2.5" resolution 
Deep image (original rms) and in the local rms map generated with a 50 pixel mesh size (corresponding to 17.5"). Contours illustrate the 
distribution of points in the most densely populated part of the scatter plot. 
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VLA-COSMOS Large (light grey curves) and Deep mosaics (dark grey curve). A cubic spline fit to the measured rms around the bump 
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panel), b) Fraction of sources in the Large catalog (vl.O) at a given radial distance r that have a S/N smaller than 4cri3(,op in the Deep 
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Fig. 7. — Fraction of sources in the VLA-COSMOS Large catalog (vl.O) which could (solid black line) or could not (dashed black line) 
be successfully assigned a counterpart within 1.5" in the VLA-COSMOS I)eep catalog. The effective success rate (grey line) is determined 
by accounting for the fraction of sources that lie below 4(t in the Deep image (cf. panel (b) of Fig. |6]) and thus cannot be included in the 
Deep catalog. The vertical axis is normalized to the total number of sources in the catalog that lie in a given concentric ring of width 2'. 
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Fig. 12. — Comparison of two diflfcrent approaches to correct for the bandwidth smearing (BWS) effect. (See text for details) 



24 



Schinnerer et al. 



10 


y 1 






1 


1 1 


8 


\ 








- 


6 


\ 








- 


4 


ji'. ... ■ . ■' 








- 


2 

1 
0.8 






••".V^ "."".■ 


.• . ." •■•■ ■■.'/ . :■ ; 


~ 


--f.,^ ----- 








- 


0.6 


y^ 








- 


0.4 


- ^ /^ 






1 





10 



100 



S/N 



Fig. 13. — Ratio of total Stotal E'ld peak flux density Speak, corr. (BWS corrected) as a function of the detection S/N, defined as the 
ratio of the (uncorrected) peak flux density and the local (17.5" box) rms. The solid lines show the upper and lower envelope of the region 
defined to contain unresolved sources (small dots), based on the simulations presented in the Appendix. Sources lying above the upper 
envelope are considered resolved (large dots). 



300 



a) 



n unresolved sources 
■ resolved sources 



200 



100 




0.05 0.1 0.5 

Sp..k, corr. [mJy/beam] 



300 



200 



100 




Fig. 14. — Number distribution of unresolved and resolved sources in the VLA-COSMOS Joint catalog as a function of corrected peak 
(top) and integrated flux density (bottom). Sources with 5peak, corr. or 5total > 1 mJy are added to the right-most bin of the histograms. 
In panel (b) the dashed line (normalized to the counts at 0.1 mjy) illustrates the decline of the integral sources counts ~S^^^^^ which is 
expected from the flat part of the Euclidan source counts present at low flux levels fe.e. IBondi et ani2008IV 



VLA-COSMOS Survey IV. 



25 



o 

S 
m 
O 
u 
I 



1 - 



0.5 



0.1 - 



0.05 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1.1 


• src. unresolved in both Large cat. v2.0 & in Joint cat. ^' /^ 


• / 


7 


.V 


• .-/ 


. -f 


.■ •*/ 


. '.**' 


'.'/ 


. • «^ 


• • • / 


. •;•':•/■ 




' Vf' ' 


. -..V. 


..' ,"^ 


:.';/;• 


• '.r.'-^ 


•r::^'^ 


» •.'••>/• 


.•*."C^ 


• 'O^?^' 


.•?V^;, 


" • *• ^i)ri^* 


. I"* ,"«v!sr; 






• ''• 'tt3^^^ * 20 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 










, . vi^;5f^^- 


! / '- 




. »j^^" '' 


'. / I 




A'^^-'. . 10 


L /J 




^mS39^^9^ 


/ 




» jMBt^jK?* " * 


•/ 


" 


' iBHEr^> * 5 


'— ■■/ — 


" 








/ W^' 


-f , , 1 , 1 , 1 1 1 1 , 1 1 1 1 , 1 , 1 r 


- 


1 1 1 1 1 1 1 1 1 1 1 1_ 


0,5 10 15 20 

III 1 1 


1 



0.05 



0.1 



peak, corr. 



0.5 



(VLA-COSMOS Large cat., v2.0) [mJy/beam] 
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Fig. 16. — Comparison of the integrated flux Stotal in the VLA-COSMOS Large (v2.0) and Joint catalog (the layout of the figure matches 
that of ll5|l . The majority of the multi-component sources (dark grey crosses) have larger measured integrated fluxes in the Joint catalog 
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Fig. 17. — Comparison of representative radio surveys conducted at 20 cm with the VLA and ATCA. Filled squares denote surveys that 
were conducted in the VLA A configuration reaching resolutions of Ri2", while surveys having resolutions of 5-10" (from the VLA B array 
and ATCA) are shown as filled circles. The VLA-COSMOS surveys are marked by open squares (Deep: top, Large: bottom), showing that 
they cover the largest area at their resolution and sensitivity. In order to reasonably compare the different surveys, we used the sensitivity 
that is achieved for at least 80% of the area covered. Thus in the case of surveys that have a non-uniform rms pattern the sensitivity value 
used here can be significantly higher than the best value present in the deepest part of the images (for detailed numbers see Tab. I?}. 
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Fig. 18. — Intrinsic source sizes of simulated Gaussian sources vs. their detection S/N. The dashed hne marks the upper expectation 
value 9obs. of the observed size of a point source of varying S/N in the Deep mosaic. At low S/N, point sources may be significantly larger 
than the beam size, 8^. The intrinsic (deconvolved) source size, flumiti inferred from the dashed line is shown in red. It defines an upper 
limit to the intrinsic source size below which a conclusive classificatio n as resolved or unresolved is no longer possible. Potentially resolvable 
simulated sources are plotted in grey for easier identification in Fig. 1191 
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In the upper panel we show the fraction of sources in a g iven bin of S/N that (1) have an intrinsic source size S < 9ii„jit but which 
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Fig. 20. — Comparison of the recovered total flux density Stout of simulated Gaussian sources with the input value, 5tin. Top row: 
distribution of relative errors (right) and plot of the relative errors as a function of sign al-to -noise ratio (SNR; left) for all simulated sources 
after these have been classified as resolved or unresolved using cq. I IA2I I. (Sec also Fig. 1191 ) Bottom row: close-up of the simulated sources 
in the bin vfith the smallest simulated S/N. 

Total flux densities tend to be overestimated by an amount that increases towards lower S/N. The according average bias is listed in Tab. 
[5] for resolved sources as well as for the complete sample of resolved and unresolved simulated sources. 
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TABLE 1 

VLA Pointing Centers for 

VLA-COSMOS Deep Project 

Pointing # R.A. {J2000) DEC (J2000) 



F07 


10:00:58.62 


+02:25:20.42 


F08 


09:59:58.58 


+02:25:20.42 


Fll 


10:01:28.64 


+02:12:21.00 


F12=' 


10:00:28.60 


+02:12:21.00 


F13 


09:59:28.56 


+02:12:21.00 


F16 


10:00:58.62 


+01:59:21.58 


F17 


09:59:58.58 


+01:59:21.58 



Note. — The naming convention 

of the VLA-COSM OS Large project 
l lSchinnerer et aT]|2007l ) has been kept for the 
pointing centers of the VLA-COSMOS Deep 
project at 1.4 GHz. 
■^ COSMOS field center 



TABLE 2 

Source Numbers in the VLA-COSMOS catalogs 



Catalog 


Type 


# of objects 


Remarks 


Deep Project 
Deep Project 


components'* 
sources'' 


3744 
3441 


4(T (~45/iJy); total 
4cr (~45/iJy); total 


Large Project (vl.O) 


sources 


3643 
1226 


4.5(t; total 
<5a 


Large Project (v2.0) 


sources 


2417 

1611 

806 

78 


5o- (~60/iJy); total 
unresolved sources 
resolved sources 
multi-component sources 



Joint 



2865 total; combined Large (v2.0) and Deep Project 

2159 detected in both Largo (vl.O) and Deep image 

392 from Large Project (v2.0) catalog 

314 from Deep Project source catalog 

2329 unresolved sources 

405 resolved sources 

131 multi-component sources (inch 43 new ones) 



Note. — See text for details. 

Direct product of running the AIPS task SAD on the Deep image 

Cleaned for multi-component and 'twin' sources 



TABLE 3 
1.4 GHz Joint Source Catalog of the VLA-COSMOS Project (abridged) 



Name 



former Name 
(in Large Project catalog, v2) 



R.A. 
[deg], {J2000.0) 



[deg] 



Dec. 
, (J2000.0) 



R.A. 
(J2000.0) 



Dec. 
(J2000.0) 



O-R.A. 



CDc 



o 

CO 



COSMOSVLADP.J09582L65+024628.1 


COSMOSVLA.J095821. 65+024628.1 


149.5902208 


2.7744972 


09 58 21.653 


+02 46 28.19 


0.25 


0.25 


COSMOSVLADP_J095821. 78+024820. 6 


COSMOSVLA.J095821. 78+024820.6 


149.5907833 


2.8057333 


09 58 21.788 


+02 48 20.64 


0.35 


0.36 


COSMOSVLADP.J09582L81+014550.7 


COSMOSVLA.J095821. 82+014550.8 


149.5908875 


1.7640944 


09 58 21.813 


+01 45 50.74 


0.38 


0.41 


COSMOSVLADP.J09582L82+014724.1 


COSMOSVLA.J095821. 81+014724.2 


149.5909333 


1.7900389 


09 58 21.824 


+01 47 24.14 


0.26 


0.26 


COSMOSVLADP_J095821.94+020707.7 


COSMOSVLA.J095821. 94+020707.7 


149.5914333 


2.1188167 


09 58 21.944 


+02 07 07.74 


0.27 


0.27 


COSMOSVLADP.J095822. 11+014058. 9 


COSMOSVLA.J095822. 10+014058.7 


149.5921333 


1.6830278 


09 58 22.112 


+01 40 58.90 


0.29 


0.27 


COSMOSVLADP.J095822. 18+014524.3 


COSMOSVLA.J095822. 18+014524.3 


149.5924333 


1.7567500 


09 58 22.184 


+01 45 24.30 


0.29 


0.30 


COSMOSVLADP_J095822. 25+013512. 3 


COSMOSVLA.J999999. 99+999999.9 


149.5927083 


1.5867500 


09 58 22.250 


+01 35 12.30 


0.44 


0.38 


COSMOSVLADP.J095822.30+024721.3 


COSMOSVLA.J095822. 30+024721.3 


149.5929250 


2.7892583 


09 58 22.302 


+02 47 21.33 


-99.00 


-99.00 


COSMOSVLADP.J095822. 57+020239.1 


COSMOSVLA.J095822. 57+020239.1 


149.5940667 


2.0441972 


09 58 22.576 


+02 02 39.11 


0.27 


0.29 


COSMOSVLADP_J095822.81+023604.3 


COSMOSVLA.J095822. 81+023604.5 


149.5950583 


2.6012000 


09 58 22.814 


+02 36 04.32 


0.48 


0.43 


COSMOSVLADP.J095822.93+022619.8 


COSMOSVLA.J095822. 93+022619.8 


149.5955667 


2.4388333 


09 58 22.936 


+02 26 19.80 


-99.00 


-99.00 


COSMOSVLADP_J095823.25+020859.4 


COSMOSVLA.J095823. 25+020859.4 


149.5969083 


2.1498417 


09 58 23.258 


+02 08 59.43 


-99.00 


-99.00 


COSMOSVLADP_J095823.27+021455.9 


COSMOSVLA.J095823. 27+021455. 5 


149.5969917 


2.2488639 


09 58 23.278 


+02 14 55.91 


0.39 


0.35 


COSMOSVLADP.J095823.67+021201.4 


COSMOSVLA.J095823. 68+021201.4 


149.5986333 


2.2004000 


09 58 23.672 


+02 12 01.44 


0.29 


0.31 


COSMOSVLADP_J095824.02+024916.0 


COSMOSVLA_J095824.02+024916.0 


149.6000833 


2.8211167 


09 58 24.020 


+02 49 16.02 


-99.00 


-99.00 


COSMOSVLADP.J095824.02+025029.5 


COSMOSVLA.J095824.02+025029.3 


149.6000875 


2.8415333 


09 58 24.021 


+02 50 29.52 


0.26 


0.26 


COSMOSVLADP.J095824.13+013836.6 


COSMOSVLA.J095824. 14+013836.6 


149.6005542 


1.6435250 


09 58 24.133 


+01 38 36.69 


0.29 


0.31 



Note. — VLA-COSMOS Joint Catalog of radio sources at 1.4 GHz representing all reliable radio sources from the Large catalog and a search for sources in the Deep mosaic 
(see text for details). All measurements were made on the 2.5" X2" Deep image. The complete tabic is available via the link to a machine-readable version above and/or at the 
COSMOS archive at IPAC/IRSA". A portion is shown here as an example of its form and content. 

" |http : //www . irsa . ipac . edu/data/COSMQS/tables/vla/ 1 
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TABLE 3 
cont. 



^pcak 


^pcak.corr. 


^total 


mis 


^M, dec 


''m, dec 


PAdcc 


Flags 








[mjy/beam] 


[mJy/beam] 


[mjy] 


[mJy/beam] 


["] 


["] 


[°] 


resb 


mult': 


memb'* 


det<= 


5.218±0.029 


5.218 


5.218±0.029 


0.029 


0.00 


0.00 


0.00 


-1 











0.201±0.035 


0.201 


0.201±0.035 


0.035 


0.00 


0.00 


0.00 








-1 





0.079±0.018 


0.081 


0.081±0.018 


0.018 


0.00 


0.00 


0.00 











-1 


0.267±0.014 


0.274 


0.328±0.034 


0.014 


1.43 


0.97 


60.80 


1 











0.179±0.016 


0.182 


0.182±0.016 


0.016 


0.00 


0.00 


0.00 














0.154±0.016 


0.154 


0.154±0.016 


0.016 


0.00 


0.00 


0.00 














0.181±0.017 


0.185 


0.341±0.057 


0.017 


3.04 


1.99 


94.80 


1 











0.135±0.026 


0.135 


0.135±0.026 


0.026 


0.00 


0.00 


0.00 


-2 





1 


1 


-99.000 


-99.000 


30.120 


0.083 


35.97 


9.90 


-99.00 


1 


1 





-99 


0.153±0.017 


0.157 


0.157±0.017 


0.017 


0.00 


0.00 


0.00 














0.093±0.019 


0.096 


0.096±0.019 


0.019 


0.00 


0.00 


0.00 











-1 


-99.000 


-99.000 


116.500 


0.033 


74.99 


12.92 


-99.00 


1 


1 





-99 


-99.000 


-99.000 


5.025 


0.026 


9.90 


3.75 


-99.00 


1 


1 





-99 


0.077±0.016 


0.079 


0.079±0.016 


0.016 


0.00 


0.00 


0.00 











-1 


0.162±0.021 


0.169 


0.162±0.021 


0.021 


0.00 


0.00 


0.00 


-1 











-99.000 


-99.000 


56.620 


0.067 


56.98 


13.93 


-99.00 


1 


1 





-99 


0.461±0.029 


0.461 


0.461±0.029 


0.029 


0.00 


0.00 


0.00 


-1 











0.141±0.018 


0.141 


0.141±0.018 


0.018 


0.00 


0.00 


0.00 















Peak flux value corrected for the bandwidth smearing effect. 

Flag for resolved sources (based on Fig. |13| l: (-2) if unresolved only in VLA-COSMOS Deep image; (-1) if resolved only in 
VLA-COSMOS Large image; (0) if unresolved in both VLA-COSMOS Large and Deep images; (1) if resolved in both VLA- 
COSMOS Large and Deep images; (2) if resolved only in VLA-COSMOS Deep image. Note that the flag value of '-2' is only 
assigned to sources that were not in the VLA-COSMOS Large Project catalog. The flag '-1', on the other hand, is used for 
sources that were classified as resolved in the Large Project catalog but are listed as unresolved in the Joint catalog. 
^ Flag for distinction between sources consisting of multiple components (1) or a single component (0). 

Flag identifying whether the source was only detected in the VLA-COSMOS Large image (-1), in both the VLA-COSMOS 
Large and Deep image (0), or only in the VLA-COSMOS Deep image (1). 

^ Flag specifying if the source was detected with S/Nu 5// > 5 only at 1.5' resolution (-1); with S/N-^^ 5// > 5 at a resolution 
of both 1.5" and 2.5" (0); or only at a resolution of 2.5", either in both the large and small scale rnis map (1), or only in the 
small scale rms map (2). 



TABLE 4 
Selected Radio Surveys at 1.4 GHz 



Field 



Area^ rms Instrument 

[deg-^] [/ijy/bcam] 



Reference 



VLA-COSMOS 


(Deep) 


1.7 


27 


VLA-A 


VLA-COSMOS 


(Large) 


1.6 


20 


VLA-A 


ECDF-S 




0.23 


8.5 


VLA-A 


SSA 13 




0.25 


16 


VLA-A 


FIR9T 




10,000 


150 


VLA-B 


FLS 




5 


23 


VLA-B 


SXDF 




1.1 


27 


VLA-A 


VVDS 




0.8 


17 


VLA-B 


ATHDFS 




0.19 


80 


ATCA 


PDS 




3.65 


100 


ATCA 



this paper 

Schinnerer et al. 2007 
Miller et al. 2008 

Fomalont et al. 2006 

Becker et al. 1995 

Condon et al. 2003 

Simpson et al. 2006 

Bondi et al. 2003 

Norris et al. 2005, Huynh et al. 

Hopkins et al. 2003 



2005 



Note. — Radio surveys displayed in Fig. 1171 The values listed were derived from figures similar to Fig. 
[3] presented in the corresponding references, except for The FIRST and FLS surveys where tlic full areas are 
used. 
^ 80% of the area covered by the survey 

Highest rms value occurring for 80% of the area covered 
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TABLE 5 

Median of relative error on recovered integrated flux densities for 
simulated gaussian sources (see appendix). 

S/N range (Siin — S'iout)/'S'iin (aH sources) (5tin — 5tout)/'5'iin {resolved sources) 

5-6 -0.05±0.25 -0.17±0.42 

6-7 -0.02±0.21 -0.28±0.39 

7-8 0.01±0.17 -0.16±0.33 

8-9 -0.01±0.17 -0.19±0.29 

9-10 -0.01±0.15 -0.20±0.23 

10-15 0.01±0.12 -0.10±0.18 

15-20 O.OlitO.lO -0.08±0.15 



Note. — Errors span the range (+/-) eontaining 2/3 of the measurements. 



